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We report improved measurements of branching fractions for charmless hadronic two-body B 
meson decays containing an u> meson in the final state. The results are based on a data sample of 78 
fb _1 collected on the T(4S) resonance by the Belle detector. We measure the branching fractions 
B(B+ -» ujK+) = (6.5±i;| ± 0.6) x 10~ B and B(B+ -» uj-k + ) = (5.7±i;| ± 0.6) x HT 6 . We give 

uK°) < 7.6 x 10~ 6 and B{B° -> lutt ) < 1.9 x 10~ 6 . 



90% confidence upper limits for B(B° 

We also obtain the partial rate asymmetries Acp = 0.06t"o'ig ± 0.01 for B ± -> LUK+ and Acp = 
0.50 



+0 f ± 0.02 for B ± 



PACS numbers: 13.25.Hw, 14.40.Nd 



INTRODUCTION 



APPARATUS AND DATA SET 



Charmless hadronic B decays play an important role 
in the understanding of CP violation in the B system. 
These decays proceed primarily through interfering b —> s 
loop penguin diagrams and b — > u tree spectator dia- 
grams. Studies of B — > ojh, where h denotes K + , ir + , K°, 
and 7T°, are important examples of such decays. Charge 
conjugates are implied unless otherwise stated. Wc also 
assume equal production of B + B~ and B°B° pairs from 
the T(4S). 



TABLE I: Measurements of branching fractions for B — > 
loK + and B + -> lotv + from CLEO, BaBar and Belle. The 
units are 10 -6 . 



Mode 


CLEO[2] 


BaBar [4J 


Belle[5] 


ujK+ 

UJH + 


3.2±f;^ ±0.8 
11.3±|;| ± 1.4 


5.5 ±0.9 ±0.5 
4.8 ±0.8 ±0.4 


9.2±|;| ± 1.0 
±.2±\% ±0.5 



Table [H lists the branching fractions from previous 
measurements 0,0, 0,0,0, which indicate some discrep- 
ancies for B + — > uK + . Naive factorization and QCD fac- 



torization approaches [fJLl yield values of B(B + — > lutt' 
consistent with the experimental results. However, these 
approaches predict B(B + — * lott + ) to be a factor of two 
larger than B(B + — * ojK + ), which is not supported by 
Belle's previous experimental results that were based on a 
29.4 fb _1 data sample 0- In this paper, we update our 
previous measurements on tuK + and loh + with a 78.1 
fb _1 data sample. We also report measurements of ujK° 
and u>ir° decay modes. 



The data sample used was collected with the Belle de- 
tector at the KEKB asymmetric energy e + e~ collider ||, 
which collides 8.0 GeV e~ and 3.5 GeV e + beams at 
a small crossing angle (±11 mrad). The data sample 
contained 85.0 x 10 6 BB pairs produced at the T(45) 
resonance. A 8.8 fb _1 data sample taken at a center-of- 
mass energy 60 MeV below the T(AS) is used to charac- 
terize continuum background. In order to establish the 
event selection criteria, we use a Monte Carlo (MC) gen- 
erator [Pj to generate signal, generic b — * c, and other 
charmless rare B decays. The GEANT3 package 01 is 
used for detector simulation. 

The Belle detector measures charged particles and pho- 
tons with high efficiency and precision llj. Charged par- 
ticle tracking is provided by a silicon vertex detector 
(SVD) and a central drift chamber (CDC) that surround 
the interaction region. The charged particle acceptance 
covers the laboratory polar angle region between 9 = 17° 
and 150°. Charged hadrons are distinguished by combin- 
ing the responses from an array of silica aerogel Cercnkov 
counters (ACC), a barrel-like array of 128 time-of-ffight 
scintillation counters (TOF), and dE/dx measurements 
in the CDC. The combined response provides K/ir sep- 
aration of at least 2.5cr for laboratory momentum up to 
3.5 GeV/c. Electromagnetic showers are detected in an 
array of 8736 CsI(Tl) crystals (ECL) located inside the 
magnetic volume, which covers the same solid angle as 
the charged particle tracking system. The magnet return 
yoke consists of alternating layers of resistive plate coun- 
ters and 4.7 cm thick steel plates for detecting if°'s and 
identifying muons. 
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EVENT SELECTION 

Hadronic events are selected using criteria based on 
the charged track multiplicity and total visible energy 
sum; the efficiency is greater than 99% for generic BB 
events . All primary charged tracks must satisfy qual- 
ity requirements based on their impact parameters rel- 
ative to the run-dependent interaction point (IP). The 
deviation from the IP position is required to be within 
±1.5 cm in the transverse direction and ±2 cm in the 
longitudinal direction. Charged particle identification is 
based on the ratio KID= Ck/{Ck + Ck), where Ck and 
C„ are likelihoods for K and tt hypotheses. A higher 
value of KID indicates a more kaon-like particle. 7r° me- 
son candidates are reconstructed from pairs of photons, 
each consisting of energy clusters greater than 50 MeV, 
with 77 invariant mass within 3er {a = 5.4 MeV/c 2 ) of 
the 7T° mass. Kg meson candidates are reconstructed 
using pairs of oppositely charged particles that have an 
invariant mass in the range 480 MeV/c 2 < to(7t + 7t~) < 
516 MeV/c 2 . The vertex of the Kg candidate is re- 
quired to be well reconstructed and displaced from the 
interaction point, and the Kg momentum direction must 
be consistent with the K s flight direction. Candidate 
lo — ► 7T + 7r~7r decays are reconstructed from charged pi- 
ons with KID < 0.9 and tt°s with center-of-mass frame 
momentum greater than 0.35 GeV/c. The uj meson can- 
didates are required to have an invariant mass within ±30 
MeV/c 2 of the nominal value (±2c). 

B RECONSTRUCTION 

B meson candidates are formed by combining an to 
meson with either a kaon (K + , K°) or a pion (ir + , 7r°). 
We require KID > 0.6 and KID < 0.4 for K+ and tt+, 
respectively. Studies from D*+ -> D Q ir + (D° -> K--K+) 
decays give particle identification efficiencies, €k = 85% 
and = 89% with misidentification rates, = 8% and 
Jk = 11%, respectively. 

B meson candidates are then identified using the beam 

constrained mass M bc = tJ{E™J 2 - \P§ M \ 2 and the 

energy difference AE = E% M - E™ m , where E^ m = 
5.29 GeV, and P§ M , Eg M are the momentum and energy 
of the B candidate in the T(4S I ) rest frame. For the AE 
calculation, the kaon in candidate B + — » uK + decays 
is assigned a pion mass so that ujK + and lot: + can be 
fit simultaneously. For events with multiple candidates, 
the best candidate is selected using the quality of the 
B vertex fit. According to signal MC, the resolutions 
for M bc and AE are 3 MeV/c 2 and 24 MeV respectively 
for B — > ujK + , ujtt + and uK° decays. For the decay 
B — > W7r°, the resolutions are 3.5 MeV/c 2 for M bc and 
55 MeV for AE. 

The B candidates are required to be within the rectan- 



gular region in the M hc - AE plane, 5.2 GeV/c 2 < M bc < 
5.3 GeV/c 2 and \AE\ < 0.25 GeV. Signal regions of 
5.27 GeV/c 2 < M hc < 5.3 GeV/c 2 and \AE\ < 0.10 GeV 
are used to display fit projections. Sideband regions 
are defined as 5.2 GeV/c 2 < M bc < 5.26 GeV/c 2 with 
\AE\ < 0.25 GeV for AE, and 5.2 GeV/c 2 < M hc with 
0.10 GeV < \AE\ < 0.25 GeV for M hc . 

Since B^w/iisaP^VP decay, where V means 
vector and P means pseudo-scalar particles, the u> me- 
son is polarized. The u> helicity angle, 0hei, is defined as 
the angle between the B flight direction and the vector 
perpendicular to the u> decay plane in the u rest frame. 
Further background suppression is achieved using the u> 
helicity and the quality of the B vertex fit (xb)- 

BACKGROUND SUPPRESSION 

Backgrounds from b — > c decays and the feed-across 
from other charmless rare B decays are found to be neg- 
ligible using MC simulations that assume the best known 
branching fraction for each decay. The dominant back- 
grounds arise from the e + e~ — > qq (q = u, d, s or c) 
continuum process, which has a jet-like event topology 
in contrast to the spherical BB events. 

Several event-shape variables arc used to distinguish 
between B decays and continuum background. The 
thrust angle 0t is defined as the angle between the pri- 
mary B decay daughter u> and the thrust axis formed by 
all the particles from the other B. S± is the scalar sum 
of the transverse momenta of all particles outside a 45° 
cone around the primary B decay daughter direction di- 
vided by the scalar sum of their momenta. In addition to 
these, a set of variables derived from Fox- Wolfram mo- 
ments are used. The moments are defined by 

RSO _ T,i, k \Pi\\Pk\Pi(cos8 ik ) 

and 

ROO = Yhj Nbj|fHcos%) 

Yh,j \Pi\\Pj\ 

where p stands for particle momentum, and Pi is the 
I th Legendre polynomial. There are two groups of par- 
ticles that go into this summation. The index k refers 
to (neutral or charged) particles from the B candidate, 
while i and j refer to other particles not from that B can- 
didate. Rf°, H3 and R°° are not used because of their 
strong correlation with M bc . To optimize the discrimina- 
tion, the remaining 5 variables (I < 4) are combined with 
cos 0t and S± to form a Fisher discriminant T 0, ^| . 
The cosine of the angle between the B flight direction 
and the beam axis (cos#b), and T are found to be inde- 
pendent, and their probability density functions (PDFs) 
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are obtained by using MC samples for signal, and off- 
resonance data for continuum background. The variables 
cos 9b and T are then combined to form a likelihood ratio 
CR = C s /(C S +Cb g ), where £~ s (bg) is the product of signal 
(qq) PDFs. A selection requirement is imposed on CR to 
reject continuum background. A typical cut is CR > 0.5 
and retains approximately 83% of the signal candidates 
while reducing the background by approximately 73%. 



ANALYSIS 

Signal yields are obtained using M\, c and AE as in- 
dependent variables in an extended unbinned maximum 
likelihood (ML) fit after restrictions are imposed on the 
variables x%i £R an d cosflhd- These are: CR > 0.65 
for uK + /tt+, CR > 0.5 for uK°, CR > 0.8 for wtt and 
| cosflhci I > 0-5- For N input candidates, the likelihood 
is defined as 



C(N s ,N B ) = e 



NB) Uli[NsP Sl (M hc )Ps^E) 
N B P Bi (M hc )P Bi (AE)}, 



where the index i runs over each event, and P Bi are 
the probability densities as functions of Mb c and AE for 
signal and background, respectively. This method treats 
the extracted yields for signal N$ and background N B 
according to Poisson statistics and constrains their sum 
to the observed number of candidates N at the maximum 
likelihood. 

The signal PDFs are determined from signal MC while 
the continuum background PDFs are derived from the 
off-resonance data. The background shapes are verified 
using data from the sideband region. The PDF for the 
AE background is modeled by a second-order polynomial 
function. The PDF for the Mb c background distribution 
is modeled with a smooth function with parameters de- 
termined from off-resonance data 16]. To model the low 
energy tail, the AE signal PDFs use a "Crystal Ball" 
line shape function |l7| with parameters determined by 
fits to signal MC. The M^ c PDFs are the sum of two 
Gaussian functions with different widths, which were ob- 
tained by fits to signal MC. Studies of B + — > D°n + and 
D° — > K + ir~ir° decays were used to fix the mean M^ c . 
Differences between widths obtained in these studies and 
those from the signal MC are regarded as systematic un- 
certainties. 

The overall reconstruction efficiencies, e, are the prod- 
ucts of detection efficiencies, determined from MC with 
no KID requirements, and KID efficiencies determined 
from D*+ -> D°7r+, D° -> K~Tr + events in the 
data. The statistical significance (S) is defined as 
y/— 21n[£(0)/£ max ], where £ max ^ s ^ e maximum likeli- 
hood at the nominal signal yield and C(0) is the like- 
lihood with the signal fixed at zero. The 90% confi- 
dence level upper limit is calculated from the equation 
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FIG. 1: Signal region projections of Mhc (left) and AE (right) 
for luK + , uj-k + , uiK° and urn . The solid curves show the 
results of the 2D fits with the background components rep- 
resented as dashed curves. Small background enhancements 
near MeV in (b) and —50 MeV in (d) are from misidentificd 
B + — > uj-K + and B + — > loK + decays. 



90%, where only the statistical uncertain- 



J^max C ( x -j dx 

ties are considered. For the final upper limit, the above 
limit is increased by one standard deviation of the sys- 
tematic error. 



TABLE II: Signal yields^s), emciencies(etot) including sec- 
ondary decay branching fractions, fit significances(E), branch- 
ing fractions(B), and 90% confidence level upper limits (UL) 
on the branching fractions for ujK° and con . 



Mode 


N s 


e to t(%) 


E S(xl0~ 6 ) 


UL(xlO" 6 ) 


LOK+ 


U.6+H 


8.1 


7.8a ± 0.6 




CJTV + 


42.1^°3 1 


8.7 


6.0cr 5.7t\i ± 0.6 




loK° 


n.ili-2 


3.3 


3.2<r 4.0±i.g±0.5 


7.6 


L0TV° 




5.2 




1.9 



MEASUREMENTS OF BRANCHING 
FRACTIONS 

The results from the fits are shown in Table |n| Fig- 
ureHshows the Mbc and AE distributions, where events 
in the Mb c (AE) plots are required to be in the AE (Mb c ) 
signal region after all selection criteria. The signal yields 
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FIG. 2: a): AE distributions without KID requirements in 
the Mbc signal region. Dotted(dashed) curves indicate the 
signal components u)K + (u>tv + ) obtained from the AE fit to 
the B + — > ujh + candidate events. The solid curve shows the 
sum of signal and continuum background components, b): 
KID vs AE distributions in the Mbc signal region. Arrows 
show the signal AE region with KID > 0.6 for coK + and 
KID < 0.4 for W7r + . 



from the fits are N uK + = 44.6+ 91 

F5.2 



8.3' Nunl 



and N uK o = 11.1_ 4 ; 4 (statistical errors only). No sig- 
nal is observed for B° — > urn . For B + — > uin + and 
uiK + , the tt + / K + feed- across is not negligible and its 
level is fixed in the ML fit. For B + — > W7r + , the contri- 
bution is estimated by using uiK + yields from the fitted 
B + — » uiK + candidate events assuming no feed-across 
from W7r + , dividing by the kaon efficiency and multiply- 
ing by the kaon mis-identification probability. The result 
is 4.8±1.0 events from B + — > ujK + in the B + — > uitt + sig- 
nal. This value is consistent with the level determined by 
repeating the fit for B + — > ujtt + with the level of u>K + 
feed-across left as a free parameter: 11.6 ± 9.0 events. 
This difference is assigned to the systematic error of ujK + 
feed-across for ujtt + decay. A similar procedure is used 
to determine the u>ir + contamination in the B + — > toK + 
yield. Here the feed-across is found to be 3.3±0.7 events. 
The final measurements of branching fractions are listed 
in Table El 

Because of the assignment of the pion mass to the kaon, 
the B + — > ujK + signal peaks at AE = —50 MeV, which 
provides some discrimination from B + — -> uin + events, 
which peak at AE — MeV. We use this to provide 
a consistency check of the uK + and luit + yields by fit- 
ting to the AE distribution for B + — > uh + candidates 
with no KID requirements applied. The signal yields 
are 60.0±^;f and 47.7±J|;f events for uK+ and wtt + , re- 
spectively, which are consistent with the results using the 
KID, where the efficiency-corrected yields are 52.4^g g 7 
events for usK + and 47.3^^'q events for uj7r + . Figure |3 
shows the results of this fit and a lego plot of KID ver- 
sus AE. From the lego plot, there is a clear separation 
in the KID distribution between luK + and lutt + signal 
yields, which also provides a consistency check. 

We also examine the properties of u> candidates in our 




7t + 7i"7i inv. mass (GeV/c ) 



FIG. 3: Fitted yields in bins of (a) tt + tv~it invariant mass 
and (b) cosine of lu helicity angle for loK + and u>-k + . Solid 
curves show the distribution from signal MC normalized to 
the results from the fits. 



fit sample. The clear oj mass peak and polarized cos 6>h G i 
distribution shown in Fig. 01a) and (b) confirm our fit- 
ted signals are from uj mesons with no significant non- 
resonant 7T + 7r~7r contribution. Several other consistency 
checks have also been performed including tightening CR 
requirements, and performing 1-D ML fits to Mh c and 
AE. All studies yield consistent results. 

Systematic uncertainties for each mode are presented 
in Table IIIII with contributions arising from the back- 
ground suppression, reconstruction and the fitting func- 
tion variations. The systematic errors from the MC mod- 
eling of the Xb an d requirements are studied with 
B+ -> D°tt + and D° -> K+ir-n decays, which are 2.0% 
for CR and 3.0% for Xb- We study the systematic error 
associated with the lo polarization (cos #i lc i) requirement 
cos#hci| < 0.5 by comparing the MC distributions with 
the fitted yields distribution(Fig. EJ). We assign a 2.6% 
systematic error. The total systematic error from back- 



TABLE III: Systematic errors for uih. Feed-across means 
loK + {-k + ) for wk + (K + ) decay. The unit is in percent (%). 



Mode 



ujK^ 



uK° 



Background suppression 4.4 
Reconstruction 8.0 

Fit t\% 
Feed-across ±1.6 

N B b 1 

B(UJ -> TT + TT-TT ) 0.8 



4.4 

8.0 



4.4 
9.5 



+ 1.6 +4.4 

-1.8 -4.3 
+3.6 
-3.3 

1 1 

0.8 



0.1 



4.4 
9.4 



1 

0.8 



Sum 



9.4 10.0 11.4 10.6 



ground suppression is 4.4%. The systematic error due 
to uncertainties in the reconstruction is determined from 
detailed studies of the charged particle tracking, uj mass 
resolution, KID and ir° detection. For charged tracking 
and 7T° detection, the decay modes r\ —> 77,7r + 7r~7r° and 
77 — » 7r 7r°7r° are used. By comparing results in data and 
MC, we assign a relative error of 2.0% for charged track 
reconstruction, 3.0% for ui mass cut, 4% for 7r° detec- 
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tion. The total reconstruction systematic error ranges 
between 8.0% and 9.5% . The mean and width differ- 
ences of Mbc and AE distributions between data and 
MC from B + — » D°ir + decays are included in the sys- 
tematic errors from fitting. The systematic uncertainty 
on the branching fraction of lu 
from the PDG tables [3. 



A C p MEASUREMENTS 



We determine partial rate asymmetries defined as 



A 



CP 



N(B~ -> uh~) - N(B+ -> uh+) 
N(B- -> ujh-) + N{B+ ^ Loh+) ' 



The values of Acp were measured for the modes £? d 



wif ± and B ± 



by performing 2D Mb c — A_E fits 



to the i? + and B separately, as shown in Fig. QJ The 
number of signal events in the ujK , and urn 1 * 1 modes are 
21.0±|;^, and 10.7±H for B+ decays, and 23.6±|;|, and 
32.2~!Lja for B~ decays, respectively The corresponding- 
partial rate asymmetry values are A C p = 0.06±g;?|±0.01 
for uK* and Acp = 0.50±£^ ± 0.02 for wtt*. The sys- 
tematic error in Acp comes mainly from the reconstruc- 
tion efficiency of high momentum charged particles and 
the fitting functions. The latter is measured by varying 
the parameters of the fitting functions. The asymmetry 
in K reconstruction efficiency is studied with an inclu- 
sive charged kaon sample. 

In the confidence level calculation, we expand the in- 
terval determined solely from the statistical error by one 
standard deviation of the systematic error. The 90% con- 
fidence level interval corresponds to —0.25 < Acp < 0.41 
for B ± -> uK ± and 0.15 < A C p < 0.90 for B ± -> wtt*. 
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FIG. 4: Projections of the 2D M bc and AE for B ± -> wK* 
and B — > ujtt^ 1 decays. Solid curves show the fit results. The 
dashed curves indicate the backgrounds. 



We also search for partial rate asymmetries in B ± — > 
wK* and wtt ± . We find A C p = 0.06±g;^±0.01 for ujK ± 

and Acp = O.SOIqJo ± °- 02 for W7r± - These correspond 
to 90% confidence level intervals of —0.25 < A cp < 0.41 



wK* and 0.15 < A cp < 0.90 for B ± -> lutt- 



for B ± 

Our results indicate the possibility of non-zero A cp for 
S ± — > W7T ± with 99.2% confidence level, equivalent to 
2.4cr significance for Gaussian errors. 
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